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ABSTRACT 

The GMRT rcioiiization effort aims to map out the large scale structure of the 
Universe during the epoch of reionization (EoR) . Removal of polarized Galactic emis- 
sion is a difficult part of any 21 cm EoR program, and we present new upper lim- 
its to diffuse polarized foregrounds at 150 MHz. We find no high significance ev- 
idence of polarized emission in our observed field at mid galactic latitude (J2000 
08h26m-|-26). We find an upper limit on the 2-dimensional angular power spectrum 
of diffuse polarized foregrounds of [PCi/2'kY/'^ < 3K in frequency bins of width 
5v = \ MHz at 300 < / < 1000. The 3-dimensional power spectrum of polarized emis- 
sion, which is most directly relevant to EoR observations, is [fc'^Pp(A;)/27r^]^/^ < 2K 
at A;_L > 0.03/iMpc~^, k < 0.1ft- Mpc~^. This can be compared to the expected EoR 
signal in total intensity of [k^P{k)/2T:^]^^^ ^ lOmK. We find polarized structure is 
substantially weaker than suggested by extrapolation from higher frequency obser- 
vations, so the new low upper limits reported here reduce the anticipated impact of 
these foregrounds on EoR experiments. We discuss Faraday beam and depth depolar- 
ization models and compare predictions of these models to our data. We report on a 
new technique for polarization calibration using pulsars, as well as a new technique to 
remove broadband radio frequency interference. Our data indicate that, on the edges 
of the main beam at GMRT, polarization squint creates ~ 3% leakage of unpolarized 
power into polarized maps at zero rotation measure. Ionospheric rotation was largely 
stable during these solar minimum night time observations. 



1 INTRODUCTION 

A current frontier in observational and theoretical cosmol- 
ogy is the search for and understanding of large scale struc- 
ture during the epoch of reionization (EoR). Detection of 
reionization at z ~ 10 would allow study of the emergence of 
the first abundant luminous objects. Such observations holds 
the promise of examination of astrophysical and cosmolog- 
ical processes at epochs as early as a few hundred million 
years after the start of the expansion of the Universe. 

The WMAP satellite has measured polarization in the 
Cosmic Microwave Background (CMB) at large angular 
scales. This polarization is believed to arise fro m Thomson 
scatt e ring of the CMB pho tons near the EoR (|Nolta et al.l 
l2008l : iKomatsu et all |2008l ). The observed optical depth 
r ~ 0.089 ± 0.016 corresponds to an instantaneous reion- 
ization redshift of Zrdon = 10.8 ± 1.4. 

One way to study the reionization transition in detail 
is by imaging redshifted 21cm emission. At redshifts above 



the EoR transition the gas is neutral and is predicted to 
glow with about a 25 mK sky brightness temperature. Af- 
ter reionization is complete this glow is absent. At redshifts 
close to the transition a patchy sky is expected. Simula- 
tions l|lliev et al.ll2008l ) suggest that with existing telescopes 
a measurement near 150 MHz may allow for statistical de- 
tection of ~ 20 Mpc patchiness in the neutral hydrogen. 
This detection would pin down the reionization redshift and 
begin the process of more detailed study of the transition. 

The 21 cm EoR signal, in the 10 milliKelvin range, 
is much smaller than the expected ~ lOK patchiness due 
to Galactic synchrotron emission and extragalactic fore- 
grounds. One of the most difficult challenges of 21 cm EoR 
astronomy is the removal of these bright foregrounds. The 
proposed strategy makes use of the very different spatial and 
spectral character of the two components of sky structure. 
Since the received frequency of the 21 cm line emission en- 
codes the distance of the source, the EoR signal consists of a 
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series of source screens that are largely independent. For ex- 
ample, if one detects a patchy layer of cosmic material at 150 
MHz due to the EoR transition, another patchy layer must 
be present at 151 MHz, but these two images come from 
distinct radial shells and the bright and dim patches in the 
two images are not expected to be aligned on the sky. The 
characteristic size of the ionization structures near the epoch 
of ionized-cell overlap is approximately 10 Mpc, which 
corresponds to 5 arc minutes in angular size, or about 1 MHz 
in radial distance for the redshifted 21cm line at 2 ~ 9. In 
contrast, the foreground emission is highly correlated from 
one frequency to the next, since most radio-bright struc- 
tures emit with a very smooth synchrotron spectrum. To 
remove the foregrounds some type of frequency-differencing 
technique will be needed. 

At frequencies above 300 MHz many radio sources with 
a synchrotron spectrum have been shown to be linearly po- 
larized, and this raises a concern for the frequency differenc- 
ing scheme. In propagation through the interstellar medium 
the polarization rotates on the sky because of the Faraday 
effect. In the polarization basis of the telescope this means 
that each linear polarization component can oscillate with 
frequency. If the oscillation period falls close to the frequency 
separation used for the foreground subtraction, incomplete 
removal of Galactic synchrotron emission may leave a resid- 
ual that masks the EoR signal. 

One way to avoid the Faraday rotation problem is to 
measure total intensity, described by the Stokes I parameter, 
using a instrument that rejects the linearly polarized Stokes 
components Q and U. However, all real instruments have 
various sources of leakage between these components. To 
design these instruments and plan future EoR observations 
one needs information on the polarized sky brightness in the 
EoR band near 150 MHz. 

Recent studies at 350 MHz l|de Bruvn et al.l 120061 ) re- 
port that the polarized component has much more struc- 
ture on arc minute angular scales than the total inten- 
sity. Polarized structure amounting to several Kelvin has 
bee n reported. Scaling using a syn chrotron spectral index of 
2.6 (Ide Oliveira-Costa et al.l 120081 ). the polarized structure 
could be tens of Kelvin in the EoR band. This means that 
instrumental leakage from Q/U to I could severely limit the 
ability to remove Galactic emission. 

This paper places new constraints on polarized sky 
brightness in the EoR range of frequencies near 150 MHz. 
We find the level of polarized sky brightness is well below 
that expected by extrapolation from higher frequencies. We 
discuss possible mechanisms which would lead to this strong 
depolarization, and suggest future techniques to differenti- 
ate between them. 



2 THE GMRT EOR PROJECT 

Our group has initiated an effort to search for 21cm struc- 
tures at the epoch of reionization using the Giant Metrewave 
Radio Telescope (GMRT) in India. 

The telescope consists of 30 antennae of 45m diameter. 
14 are designated "central core" antennae and lie within a 1 
km area. The remaining 16 are along 3 arms of lengths up to 
10km, designated East, West, and South. The dense layout 
of the core allows high brightness sensitivity, which is needed 



for the search for reionization. For this experiment, the 150 
MHz feeds are used. These consist of orthogonal pairs of 
folded dipoles, backed by a ground plane. These antennas 
couple to X and Y linear polarizations, but the two signals 
pass through a hybrid coupler before entering the amplifiers. 
For each antenna this results in a pair of right and left circu- 
larly polarized signals which are amplified, upconverted and 
transmitted optically to the receiver room. Later processing 
allows measurement of the full set of Stokes parameters. 

To attempt the EoR experiment, and for other work at 
GMRT, we built a new signal processing system for the tele- 
scope, the GMRT software backend (GSB, Roy et al 2009, 
in preparation). This consists of 16 commercial Analog to 
Digital sampling boards installed in an array of off-the-shelf 
computers. The AD boards have 4 input channels, and are 
connected to a common clock and trigger signal, allowing 
synchronous sampling. The sampling system is capable of 
32 MHz bandwidth but for this experiment, only 16 MHz 
was used. The passband was defined by an IF filter. All 
60 signals were sampled with 8 bit precision at 33 MSam- 
ple/s. Each AD board transfers the digitized data into its 
individual host computer. These streams are Fourier trans- 
formed in blocks of length 4096 samples at 16 bit precision. 
The 2048 complex Fourier coefficients are then rescaled to 
4 bits, and sent over a gigabit network to correlation nodes. 
Each block of 128 frequencies is sent to one of 16 software 
correlation nodes. These products, which we call visibilities, 
are accumulated for 1/4 second in 16 distinct gates. Thus 
the initial visibilities have a frequency and time resolution of 
7.8 kHz and 1/4 second (without accounting for the gates), 
respectively. 

The "gates" are essential to our polarization calibration 
technique. These gates in time are synchronized so each cov- 
ers one of 16 segments of the pulsation cycle of a pulsar in 
the field. The pulsar is a known source of polarized emis- 
sion. By comparing pulsar-on to pulsar-ofi' we can measure 
the system gain directly using a sky source. The raw 1/4 
second averaged visibilities were stored on disk, as well as 
16 fringe-stopped gated visibilities integrated for 16 seconds, 
and averaged over frequency into 128 frequency channels. 

All these signal processing calculations occur simul- 
taneously, and are structured as individual asynchronous 
pipeline processes. The processor and network speeds are 
sufficient that each calculation is completed in real time and 
there is less than 10% data loss. 



3 DATA 

Our target field was centered on pulsar PSR B0823-f26 
(J2000: 08 26 51.38 +26 37 23.79). The field was chosen to 
contain a bright pulsar, and minimal other bright sources. 
The galactic latitude is 6 = 30 degrees, and happens to be 
a relatively cold spot on the sky. 

The brightest radio source within the primary beam, 
and our primary fiux reference is 5C 7.245, located 12 arc 
minutes from the pulsar. Its fiux is estimat ed to be 1.5 Jy 
at 15 1 MHz from the 7C Lowdec Sur vey (Waldram ct alj 
Il996i ). It is a radio galaxy at 2 = 1.61 llWUlo tt ct al. 200l|). 
The source has two components separated by 16 arc seconds, 
and is resolved by the longest baselines at GMRT. 

A bright source outside the main beam, 3C200, is 2.5 
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degrees to the north. It is attenuated by the primary beam 
profile, but shows up as the strongest source in the raw po- 
larized maps due to I to P leakage in the side lobes. 

In Table [T] we list the 10 brightest sources in the field 
after the primary beam power attenuation effect from the 7C 
Lowdec Survey. Al so listed are their fiuxes from the VLSS 
survey at 74 MHz (jCohen et al."2007l)_ and the Te xas survey 
of radio sources at 365 MHz ( Dougla s et al.lll996h . Based on 
these three surveys we calculate the spectral indices and the 
errors, as listed in Table [T] 

This region was map ped at 150 MHz 
l|Landecker fc Wielebinskil Il970l ) and found to have 
brightness temperature 170K. Ide Oliveira-Costa et al.l 
l)2008t ) have reconstructed a variety of sky maps and used 
them to estimate spectral indices. From these reconstruc- 
tions, we estimate a spectral index of 2.6 for this region, so 
our sky brightness temperature varies from 200-150K over 
our spectral band (140-156 MHz). 

Data was taken from Dec 7-Dec 18, 2007, starting at 
10PM local time, going until 7AM. This period has minimal 
radio frequency interference. The data was recorded in one 
hour "scans". The analysis in the paper is primarily based 
on data taken during the night of December 9, which was 
fully reduced through the described pipeline at the time of 
writing of this paper. 



4 INTERFERENCE REMOVAL 

Terrestrial Radio Frequency Interference (RFI) was removed 
in two stages. First, spectral-line RFI was flagged. In each 
frequency subset of 128 channels, the distribution of intensi- 
ties are calculated. The upper and lower quartile boundaries 
are used to estimate the amplitude of the noise. Any outliers 
further than 3cr are masked. Each mask is 8 kHz wide, and 
4 seconds long. Approximately 1% of the data is flagged by 
this process for deletion. 

A bigger problem at GMRT has been broadband RFI. 
This is particularly problematic for short baselines. The ra- 
tionale for the our approach to broadband RFI segregation 
is that sources on the sky produce visibilities that oscillate 
as the Earth turns. Meanwhile, sources on the ground, al- 
though their intensity typically varies with time, remain at 
flxed lag for each baseline. We used a singular value decom- 
position (SVD) to sort among the terrestrial and celestial 
classes. 

We assemble the visibilities in a two dimensional rect- 
angular matrix, Vi{u,t) which is a function of baseline i, 
frequency v and time t. We find that Earth-fixed broad- 
band sources tend to flicker synchronously in all baselines 
at all frequencies, and are therefore factorizable as Vrfi = 
Li(v)T{t). We call Li{v) the visibility template, and Ti{t) the 
temporal template. The RFI visibility can thus be written as 
a sum of individual RFI sources a, 

v^RPi(i',t) = ^ifMT^"W. (1) 

a 

Each L{h')T{t) product function is a Singular Eigenvector of 
the matrix. 

In contrast, objects on the celestial sphere exhibit fringe 
rotation at a unique rate for each baseline, and thus do not 



factor. Celestial sources produce very small eigenvalues un- 
der a Singular Value Decomposition. 

A decomposition of the form ([1} can be accomplished 
via a Singular Value Decomposition algorithm as long as 
the eigenvectors are orthogonal. A pair of visibility vectors 
due to two physically separate sources will tend to have little 
overlap, since the dimensionality of this space is huge. There 
are of order SO'^ = 900 baselines between antenna pairs, each 
with four polarization combinations. Two distinct sources 
may accidentally fall at the same lag for one pair but will 
be separated from each other when other pairings are con- 
sidered. The time templates might not be as orthogonal. In 
particular electric arcing from AC transmission lines, which 
is likely a the cause of a substantial subset of the broadband 
RFI, might occur more often at the AC waveform peak. 
This timing can be common to many sources. For this work, 
only about 100 of the largest SVD eigenvalues are flagged 
as RFI, in a matrix with about 10^ entries, so even if the 
eigenvectors are not perfectly orthogonal, the very sparse 
space spanned by the first 100 eigenmodes is still given by 
an SVD decomposition. 

To carry out the procedure we begin by organizing the 
data in a matrix. Each hour scan has 14400 time records, 
which we call the rows of our matrix, the columns num- 
ber 2048x60x61 entries which correspond to the number of 
frequency channels and baselines. Each matrix entry is a 
complex visibility. This 7495680x14400 matrix is then SVD 
decomposed. The first 100 right eigenvectors are tagged as 
"RFI", and removed. This removes 0.7% of the matrix en- 
tries. 

We have also used the SVD procedure to physically 
locate RFI sources. We selected several of the brightest 
eigenmodes for further study. Near field imaging using these 
eigenvectors has allowed us to localize the sources on the 
Earth's surface with 100 meter precision. Visiting these lo- 
cations with radio direction finding equipment we have local- 
ized candidate emitters in several cases. The two brightest 
appear to originate from extraneous thin wires hanging ver- 
tically from high tension power lines. Other sources appear 
related to intermittent radiation from transformers. So far, 
we have visited only a few sources. 

The brightest SVD eigenmodes do appear to represent 
RFI, not celestial sources, however there are some limita- 
tions to the technique. The antennas slowly turn as they 
track objects on the sky and this rotation leads to phase 
changes for terrestial RFI sources. The telescope feed moves 
by half a wavelength every 12 minutes. For an RFI source lo- 
cated between two antenna, the relative delay changes by a 
full wavelength. We believe this sometimes causes individual 
sources to break up, occupying several RFI modes. However, 
as long as the source is bright, it is still identifiable, even if 
split into components. 

Another limitation: the SVD analysis tends to mis- 
identify, as RFI, celestial sources with fringes that rotate 
slowly. A baseline that is oriented perfectly North-South 
does not produce any fringe rotation, so baselines with small 
V are more prone to erroneous projection of celestial sources 
onto the RFI modes. We find that setting the RFI cut at 100 
eigenmodes has a noticeable impact only on baselines with 
\v\ < lOA. Increasing the number of cut eigenmodes widens 
the impacted v range, while substantially decreasing number 
of cut eigenmodes leaves significant residual RFI power. 
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Source 


RA [degj 


Dec [clegj 


^7C-BEAM [Jy] 




F74 |Jy| 


F365 [Jy] 




d 


Aa 


3C 200 


126.855 


29.3128 


2.135 


13.435 


26.78 


6.407 





93 





07 


5C 7.245 


126.487 


26.7331 


1.511 


1.536 


2.95 


0.731 





89 





04 


B2 0825+24 


127.172 


24.6107 


0.992 


2.922 


10.58 


0.671 


1 


77 





08 


B2 0819+25 


125.560 


25.6421 


0.951 


1.701 


3.15 


0.845 





84 





04 


PKS 0832+26 


128.798 


26.5764 


0.706 


2.124 


3.43 


1.384 





62 





10 


B2 0829+28 


128.082 


27.8792 


0.682 


1.634 


2.64 


1.030 





63 





08 


5C 7.111 


124.785 


26.2457 


0.670 


1.786 


3.44 


1.086 





82 





20 


B2 0832+26A 


128.770 


26.1785 


0.638 


1.960 


3.48 


1.156 





75 





11 


B2 0828+27 


127.841 


26.9659 


0.625 


0.889 


1.63 


0.400 





87 





03 


B2 0816+26B 


124.819 


26.7014 


0.569 


1.417 


2.82 


0.657 





94 





05 


5C 7.223 


126.028 


26.4672 


0.515 


0.584 


1.060 


0.306 





81 





06 



Table 1. The apparent brightest 10 point sources in the observing field, listed in descending order of their primary beam attenuated 7C 
Catalog flux, listed in the ¥ic—BEAM column; F70 indicates their original flux measured from the 7C Radio Survey at 151 MHz, F74 
from the VLSS survey at 74 MHz, and Fsgs from the Texas Survey of Radio Sources at 365 MHz. a is the spectral index calculated from 
the three surveys and Ao the respective errors on the spectral indices. 



Empirically, a cut at 100 eigenmodes works very well, 
most visible broadband RFI had been removed from the 
data after this process. 



5 CALIBRATION 

Polarization calibration is particularly challenging at low 
frequencies. The primary beam is wide, telescope gain is low, 
interference is frequent, the ionosphere is variable and the 
sky offers few polarized point sources. Normally, an iterative 
procedure is used, where one starts with a guess for the sky 
model, determines the system gain and phases, makes a new 
map, and uses this as a model. Errors in the sky model are 
entangled with errors in the telescope system calibration. 

We used a novel approach by observing a field contain- 
ing a bright pulsar. By subtracting the off-phase from the 
on-phase, one can lock in to the modulated emission. Of the 
16 gates, we took the gate containing the pulsar on state, 
and subtracted the average of the two flanking gates. Af- 
ter this subtraction, the measured visibilities contain only 
the pulsar. Other celestial or terrestrial signals remain only 
if they happen to be modulated in sync with the pulsar. 
While the pulsar flux may vary from pulse to pulse, each 
antenna sees the same flux, and the pulsar serves as a com- 
mon source allowing the relative gain of the 30 antennae to 
be measured. The pulsar also serves as a guide star. Pulsar 
fluxes can be variable in time, but their positions are very 
stable. 

This pulsar is strongly polarized for individual pulses 
but the polarization angle varies from pulse to pulse. The 
polarization fraction for individual pul ses has been measured 
to be 40% (|Hankins fc RankinI liooj ) at 430 MHz. Many 
pulsars have polarization fractions that decrease with fre- 
quency so the fraction may be higher at 150 MHz. Our data 
was stored as 16 second averages, which corresponds to 30 
pulses. We found the average polarization fraction for these 
averages was 20%, with an RMS scatter of 8 degrees. The 
polarization fraction decreases as longer averages are assem- 
bled. This is consistent with a highly polarized pulse whose 
angle varies from pulse to pulse, with only a slight alignment 
preference. This random angle behavior is also observed at 
430 MHz. Figure [1] shows the 9 minute averaged angle, after 
correcting for parallactic angle rotation. 




-2 2 4 

local time (relative to midnight) 



Figure 1. Apparent pulsar polarization angle. The solid, dot- 
ted, dashed, long dashed are the mornings of the 18th, 17th, 9th 
and 16th respectively. The vertical axis is twice the polarization 
position angle, i.e. tan~^ (fZ/Q). The long term change seen in 
the data (the dip towards the center of the plot) is likely due to 
change in airmass and geomagnetic field angle, while the variation 
between each point and those immediately following is consistent 
with pulsar angle variations at the source. The night of the 18th 
showed a two hour period in which the polarization angle departed 
from the slowly changing behavior characteristic of the rest of the 
data set. This is likely an ionospheric anomaly. For this plot an- 
gles are arithmetically averaged over one minute intervals, then 
the median angle for each 9 minute interval is plotted. 



We also used the GMRT noise injection system to mea- 
sure the gain of the electronic signal chains. At each prime 
focus a noise waveform is injected into the two polarization 
channels, modulated at 0.5 Hz. We used the noise source on 
a setting which injects 40 K, common between the L and R 
channels. We measured the strength of the noise source in 
the L-R correlation of each antenna and used this as a gain 
reference. 
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The goal of the calibration process is to determine two 
complex gains for each of the 60 signals. One is the response 
to its own polarity, and the other is the leakage from the 
orthogonally polarized mode. The pulsar model appears to 
work well, with typical RMS phase closure errors around 
1 degree averaged over 20 minutes. This model reduces the 
1830 complex visibilities into 120 complex gains for each fre- 
quency channel. The visibilities are reproduced from this re- 
duced degree of freedom model. This is called the 'bandpass 
calibration'. In addition, the model allows for 1792 complex 
gains, one for each time interval of approximately 16 sec- 
onds, corresponding to 8 hours of integration. The model 
assumes that the gains factor as a function of time times 
a function of frequency. The 1792x1830 numbers at each 
frequency are reduced to 1792-1-120. Polarization calibra- 
tion was achieved by using the parallactic angle rotation 
to decompose the pulsar into three components: the "un- 
polarized" source which is invariant under parallactic angle 
rotation, and two rotating components. More details can be 
found in Pen et al (2008). 

The system is fully calibrated using only the pulsar and 
noise sources, except for an overall scale described below. 
We find that the unpolarized pulsar flux varies by about 2% 
when averaged over one hour, relative to the average of all 
noise injection sources. The noise source and system gains 
appear stable for most antennae, with a few exceptions. The 
pulsar flux can be readily subtracted from the time averaged 
data. 

This calibration procedure results in a polarization map 
referenced to the pulsar polarization angle at each frequency. 
We t hen correct this with the pulsar's known RM=5.9 rad 
m"^ (|Manchesteilll974l : [Manchester et al.ll2005^ for the rest 
of this paper. 



6 FLUX SCALE 

We use the pulsar for relative antenna calibration, but not 
to establish the overall gain, since the pulsar flux varies 
with time. To accomplish the overall flux calibration we use 
four well-measured bright continuum radio sources in the 
field. We estimate the flux of the four brightest sources that 
lie within the full-width-half-max (FWIfM) of the primary 
beam using the estimated spectral indices of these sources 
from the catalogs listed in Table[T]and correcting for the pri- 
mary beam power attenuation. We calculate the source flux 
values from our data at 151 MHz and match the average flux 
to that of the 7C fluxes. The estimated flux and respective 
errors are listed in Table (2] The resulting fractional errors 
on these four source fluxes are less than 13%, and the error 
on the flux scale in the calibration procedure is 6%. 

The calibration of the power spectrum amplitude de- 
pends on the point source flux calibration described above 
and also on the primary beam width and profile. Through- 
out, we approximate the primary beam as a round Gaussian 
beam with a FWHM of 3.3° . This profile is a good fit to our 
holographically-measured primary beam power out to about 
2 degrees radius. 
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Figure 3. Apparent polarization of sources at Faraday depth 
(/>=0 rad/m^. The small squares are the apparent polarizations. 
The asterisks are multiplied by the primary beam, which is an 
estimate of their impact to the total polarization in the image. 
We see that most of the apparent polarized flux contribution is 
from the periphery of the beam. 



7 MODELLING 

After the primary calibration, the dominant polarized source 
is the pulsar. The second brightest apparent polarized source 
is 3C200, 2.5 degrees to the north, with an intrinsic fiux of ca 
13 Jy at 150 MHz. The source is probably not significantly 
intrinsically polarized, but rather because the side lobes are 
strongly polarized, it leaks into the polarized maps. 

Figure [2] shows the apparent luminosities and polariza- 
tions of the brightest point sources in the field. Objects ap- 
pear more polarized as they are located away from the image 
center, suggesting the primary beam as the origin of this 
apparent polarization. Figure [3] shows the radial apparent 
polarization. 

Precise cleaning of the polarized image shown in Figure 
[5] will require detailed knowledge of the off-axis leakage at 
GMRT, which we have not yet measured. Despite this im- 
perfect knowledge we have been able to make some progress 
using a modified version of a standard 'clean' procedure. We 
fringe stop the data iteratively to sources in descending or- 
der of flux as predicted from the 7C catalog convolved by 
the Gaussian primary beam model. At each source position, 
we measure the equivalent point sources flux by summing all 
baselines. The point source is modelled as two complex num- 
bers, one for the intensity and one for polarization. Separate 
values are fit every 4 MHz and every 20 minutes. A point 
source with this flux is then subtracted from the data. We 
repeat this for the 31 brightest sources in the field. For the 
short baselines, sources are confused, so we downweighted 
the central square antennae by a factor of 10. A baseline 
containing one central square antenna gets lOx less weight, 
and a baseline with two such antennae is downweighted by 
lOOx. Similarly, we downweight the last two antennae on 
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Source 


RA [deg] 


Dec [deg] 


F7C [Jy] 


Fi5l [Jy] 


Error 


5C 7.245 


126.487 


26.7331 


1.536 


1.74 


0.13 


B2 0819+25 


125.560 


25.6421 


1.701 


1.51 


0.11 


B2 0828+27 


127.841 


26.9659 


0.889 


0.93 


0.05 


5C 7.223 


126.028 


26.4672 


0.584 


0.53 


0.09 



Table 2. Sources in the observing field that are used to estimate the overall calibration flux scale. Fjc is again the flux from the 7C 
catalog at 151 MHz, Fisi is the flux estimated from our data at 151 MHz, and error is the respective fractional error on the estimated 
flux. 




-2 





RA (degrees*cos((5)) 



2 



Figure 2. Polarized point source map at Faraday depth = rad/m^. The radius (mean of semi-major and semi-minor axes) of each 
ellipse is proportionate to apparent flux. EUipticity is proportionate to polarization. Map is averaged over 16 MHz. The asterisk in the 
center is the pulsar. 5C7.245 has been omitted for clarity. The ellipticity is plotted eight times larger than the apparent polarization. The 
global polarization angle was rotated to align the pattern radially. The feed support structures could lead to such a polarization pattern 
(see text). The dashed ring is the half power point of the primary beam. Sources outside that radius appear more strongly polarized. 
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each arm, since some of the sources are resolved by them. 
This effectively uses the intermediate baselines to solve the 
model, and uses the cross correlations from arm antenna to 
central square to extrapolate the flux on the shortest and 
longest baselines. 

Each of the 31 point sources is fitted with 80 parame- 
ters: 4 frequencies by 20 temporal fluxes (one per 20 min- 
utes, total of 400 minutes), for a total of 248 degrees of free- 
dom. To quantify the amount of spurious signal removal, we 
added a 30 mjy trace in the form of a Gaussian footprint 
centered at u=100, v=-100 with a RMS of a = 20 in the 
Fourier space. Applying the same CLEAN on this data, and 
subtracting the original CLEAN data recovers the original 
footprint to better than 99%, while leakage into originally 
empty baselines acquires up to 6% of the peak flux. We con- 
clude that that this process has at most a small impact on 
the flux at \u\ ~ 140, i.e. I ~ 900. 



length is removed, from a 45m dish. The effective area of 
the dish is about lOOOm^, of which 10% is blocked, and this 
blockage is different for the two polarizations. This could 
make the E and H beams on the sky different. The blockage 
attenuates the main lobe of the beam, while generating an 
extended pedestal in the pattern. In such a model, one ex- 
pects unpolarized objects within the first null of the primary 
beam to appear radially polarized (in terms of the electric 
field). Beyond the null, objects would be tangentially po- 
larized. Another source of beam squint occurs because the 
illumination of the dish by the feed is not exactly round. 

In principle, one could measure the primary beam for 
both polarizations and correct the leaked maps to produce a 
clean Stokes I map. This would require a non-local calcula- 
tion, both on the sky and in Fourier space, and would require 
a the use of full m atrix inversion. Codes like CBI-GRIDR 
l|Mvers et al.ll2003l ) could be adapted for this purpose. 



8 ROTATION MEASURE SYNTHESIS 

At low frequencies, Faraday rotation can have a large effect. 
The pulsar has a rotation measure of RM=5.9 rad/m^. The 
apparent polarization angle rotates as x = RMA^-|-constant. 
We use units of radians per meter squared throughout this 
paper. At 2 m wavelength, the angle changes 5.3 radians 
over the 16 MHz observed band. Since polarization angles 
are periodic over 180 degrees, this results in almost two po- 
larization rotation periods. After pulsar referenced polariza- 
tion calibration, all angles are defined relative to the pulsar. 
A direct 2-D map would strongly suppress sources whose 
RM is more than ARM> 3.5 away from the pulsar's, since 
the rotating polarization angle causes strong cancellation of 
the source flux. 

To make maps at various Faraday depths, we ap- 

! ily a method known as rotation measure synthesis 
Brentiens fc de Bruvnll2006h . The frequency channels are 
mapped into a 3-D grid of u — v — . We chose a 
2048x2048x128 grid, where all our frequency channels are 
scaled to flt in the flrst 64 channels. Padding with 64 
zeros increases the sampling in Faraday depth and results 
in several reconstructed slices per RMTF width (although 
note that this does not affect the width of the RMTF). A 
3-D FFT then makes a 3D map with axes [x, y, </>) The fre- 
quencies are binned to the nearest grid point, resulting in a 
2/-K = 64% bandwidth depolarization at the Nyquist limit 
of (/> = 56 rad m~^. 

At Faraday depth = 0, we can clearly see the beam 
squint: point sources in the side lobes appear radially po- 
larized. Our calibration procedure equalizes the polarized 
response at the pulsar position, but the leakage will vary 
across the beam in a quadrupolar pattern, producing the 
radial ellipses in Figure (2] This squint, a difference between 
the width of the E and H plane antenna patterns, can arise 
from a variety of reasons. For example, the prime focus in- 
struments are mounted on four support legs, which are each 
narrower than our observing wavelength (2m). They carry 
reflecting currents more efficiently along the legs than trans- 
verse. The illumination pattern in one linear polarization has 
a horizontal line blocked, while the orthogonal polarization 
pattern would have the opposite line blocked. To estimate 
the size of such a squint term we assume of order one wave- 



9 RESULTS 

We compute the angular power spectrum at different rota- 
tion measures. The visibilities were directly binned into a 
power spectrum assuming a 3.3 degree Gaussian beam. In 
this case, 

where V is the visibility in Jansky for a baseline u measured 
in wavelengths, and 9i, is the FWHM of a Gaussian beam. 
The left hand side is measured in . It follows that I = 2'ku. 

In order to minimize the impact of non-coplanar w 
terms, we only used baselines with |w| < 200. We convolved 
each visibility with a Gaussian of RMS cr = 23 in the u-v 
space to reduce the field of view to the central square de- 
gree, which further reduces the impact of non-coplanarity, 
and the beam squint. The phase errors from non-coplanar 
w terms are /^(j) ~ wA9^ /2, and we expect phase errors up 
to a few degrees at the edge of the reduced field. 

A given point in u-v space is typically covered by several 
baselines, at different frequencies. The RM synthesis maps 
are a sum of these points, weighted by the phase factor that 
accounts for the oscillation due to the chosen rotation mea- 
sure. Power spectra can be used to summarize the statis- 
tics of these maps (Appendix . There are several types 
of power spectra that can be computed, both 2-D and 3-D, 
containing different types of information. Perhaps the most 
straightforward is the square of each gridded visibility in a 
frequency interval v, 

Vlr = Q + iU is the complex L-R visibility which contains 
both Q and U parameters. This is a conservative overes- 
timate of true polarization on the sky, since this statistic 
includes thermal noise, polarization leakage, RFI, etc. By 
construction, this power is positive definite. We measured 
in frequency bins Av = 1 MHz, so Faraday rotation within 
a bin is not expected to play an important role. 

Since different baselines see the sky at different par- 
allactic angles, one can also measure the cross correlation 
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between distinct visibilities. To reduce the impact of leak- 
age, RFI and thermal noise, we can discard the self-terms in 
the power spectrum. In equation ^ we note that the visi- 
bility location it = it(6, t) is a function of the label of base- 
line b and the time of observation t. Each gridded visibility 
= b' 1 1' ''^("(^'*))^*('^(^''*')) can contain con- 
tributions from different baselines taken at different times. 
Instead of summing and squaring, we can square, and drop 
all contributions where 6 = 6'. 

Computationally it is expensive to store all baselines 
contributing to a u-v grid point, so we randomly separate 
baselines into two bins, grid the visibilities from each bin, 
and compute their cross correlation. We repeat this many 
times to include all possible pairs of baselines contributing 
to a product, in this particular case 129. This results in an 
average power {C;), as well as a standard deviation 5Ci. The 
latter is a rough estimate of the standard deviation on (Cj). 
The variance on the mean is smaller than the variance on 
each subsample. If the baselines were densely distributed, 
the variance is overestimated by a factor of 2. The actual 
baselines are not dense, so sometimes one of the elements 
of the pair will be empty. This causes the variance to be 
further overestimated. The average (Cj) appears consistent 
with zero, so we roughly interpret 5Ci as an effective 2a 
upper bound, shown in Figure 21 

If one fixes a rotation measure, all frequencies can be 
summed into a single u-v grid point, which further improves 
sensitivities, and localizes the RM of the power spectrum. 
Before subtracting the point sources, excess power is vis- 
ible at (/) = when imaging a 4 degree field. After point 
source subtraction, we saw no rotation measure at which 
there is a statistically significant signal. Restricting the field 
to one degree using a primary beam convolution in u-v space 
should further suppress any potential contamination. We 
apply 6b = 1 degree in Equation (|3]|. We show two rep- 
resentative lines in Figure |3] at </> = — 14 and 14 rad m~^, 
where we have plotted the absolute value of the power. Be- 
cause we discarded auto-correlation terms, noise can lead 
to positive or negative correlations. The effective width of 
the window in Faraday measure space is determined by the 
16 MHz bandwidth, which corresponds to bins of width 

~ 3.5 rad voT^ . 

[Note: there are several ways to measure the width of 
the RMTF. For a top-hat range of frequencies, |RMTF| (</>) = 
|sinc(<^AA^)| (the phase depends on the wavelength to which 
one de- rotates). The width relevant for power spectrum anal- 
ysis (Appendix ig is / jRMTFp d</> = 7r/A(A^) since if 
a slice is suppressed by the RMTF then its power spec- 
trum is suppressed by RMTF'^. One could alternatively com- 
pute the FWHM of |RMTF|, which for the sinc-function is 
3.79/A(A2).] 

The observed apparent polarized power could either be 
polarized structures on the sky, or residual leakage from 
beam squint. Towards the edge of the beam, point sources 
appear in the polarized maps at up to 10% of their unpolar- 
ized fiux, with a quadrupolar pattern characteristic of beam 
squint: objects on two diametric sides of the fields appear as 
positive sources, and at 90 degrees appear negative. 

We show a map of the point-source-subtracted sky in 
Stokes Q and U in Figure [S] at — 3radm~^, half way to 
the pulsar RM. Excess power is still discernible in the region 
of the primary beam. The point source subtraction reduces 



the impact of this leakage, but it could still leave residuals. 
Our apparent measured polarized signal might not be truly 
polarized on the sky, but rather leakage from the unpolarized 
sky brightness or residual man-made interference. We treat 
the measured power in this image as an upper limit to the 
polarized sky brightness. 

9.1 3D Power Spectrum 

EoR structure is expected to be three-dimensional (3D), so 
we need to compare the 3D power spectrum of the po- 
larized foreground to that of the expected signal. Theo- 
rists have traditionally written the EoR power spectrum in 
te rms of the temperature variance per logarithm ic range in 
k (iBharadwai fc Alill2004l : iBarkana fc Loebll2005l ): 

Ai.(fc,M) = |^PT(fc,M), (4) 

where ^ is the cosine of the angle between the wave- vector 
and the line of sight and Pt is the 3D temperature power 
spectrum. It is also common to work in terms of the radial 
and transverse components of the wavenumber, fc|| = kfj, and 
k± — k\J\ — /i^. 

It is appropriate to evaluate EoR foregrounds by pro- 
jecting the Galactic polarization structure onto the cos- 
mological volume using the usual mapping v z = 
1420MHz/i^ — 1, and estimating the polarized 3D power 
spectrum that one would infer. Details of this power spec- 
trum definition are given in Appendix [X] At our central 
frequency v = 148 M Hz, and with the W MAP determi- 
nation of rim = 0.26 (|Dunklev et al.l |2008| ) the constants 
described in the appendix are: the redshift z — 8.59; co- 
moving radial distance D = 6720/i~^ Mpc (where Ho ~ 
100/ikms-^Mpc-i); H/c = 5.04 x IQ-^/iMpc-^ and C = 
12.9/i~^ MpcMHz"^ The mapping of Eq. (fX38|l between 
the Faraday depth and the apparent radial wavenumber is 
fc|| = 0.0086/1 Mpc"^ (radm~2)~^ (f>. The 3D polarized power 
spectrum Pp{k,fi) is simply related to the 2D power spec- 
trum Cf of the Faraday depth slice in the RM synthesis 
map at this value of 0. One can then obtain 

^Uk) = ^Ppik)- (5) 

The 3D power spectra Ap(fe) are shown in Figure 6; 
the power per In k is at most ~ 2 K at A: ^ O.l/i Mpc"'^. This 
rotation measure corresponds to polarized emission at a line- 
of-sight wavenumber of fc|| — 0.12/i Mpc"'^. The implied 3D 
power spectra are shown in Figure [S] 



10 DISCUSSION 

If extended polarized emission is bright at 150 MHz, leak- 
age from Stokes Q ox U into Stokes / may make removal 
of this foreground from 21 cm EoR images quite difficult. 
Removal of the synchrotron intensity foreground from EoR 
images requires projecting out structure that varies slowly 
with frequency (i.e. has fc|| « 0). This also removes any com- 
ponent of the EoR signal that varies slowly along the line 
of sight, i.e. it suppresses modes with |A;||| less than a few 
times {ClS.v)~^ . The anticipated EoR structure is at scales 
of order k ~ 0.1/iMpc~^; this range of scales that will be 
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Figure 4. Power spectrum of unpolarized (l'^ , long dash) and total polarized emission {P^, dotted) over the central square degree, 
after subtracting the 31 brightest sources in the primary beam. This polarized signal includes leakage and thermal noise. The upper 
bound symbols are a 2-a estimate of polarization bound using baseline cross correlation. They are averages of the polarization power in 
1 MHz bands. The lower solid and dot-dashed lines are the total power polarization measured at </<=-14 and 14 over a width A0=3.5, as 
determined from the map power. Using the cross correlation technique described in the text, we plot a less conservative estimate with 
circles for = —14 and squares for (j> = +14. Open symbols denote negative power, and filled are positive values. For details on the unit 
conversions, see the appendix. 



contaminated hy Q,U ^ I leakage of polarized structure 
at Faraday depth of order ~ lOradm"^ [see Eq. (|A38|l ]. 
This is the Faraday depth range for which polarization leak- 
age causes the most difficulty for an EoR search, although if 
the polarization leakage varies rapidly with frequency other 
values of are also of concern (see Appendix I A4p . 

Previous work at 350 MHz (|de Bruvn et al.ll2006l ) in- 



dicated polarized emission with spatial variations of several 
Kelvin on angular scales of tens of arcminutes. This polar- 
ized structure is reported to be brighter than the unpolarized 
component at the same angular scale. Extrapolating this po- 
larized structure to 150 MHz with a spectral index of 2.6, 
this corresponds to several tens of Kelvin in the EoR band. 
To remove this foreground would require that instrumental 
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Figure 5. Polarized sky image at </> = 3radm~^. The left panel is the Stokes Q intensity, while the right is Stokes U. All baselines with 
\w\ < 200 are included. The field of view is 9.78 degrees, which is three primary beams. Real polarized sky structure as well as leakage 
I into Q,U will have higher contrast within the main beam region. Sources that generally have the same contrast inside and outside the 
main beam area include RFI and the aggregate extended image artifacts of multiple bright sources far outside the main beam. These 
images do show a weak enhancement of contrast in the main beam region, but the enhancement is too weak to allow confidence in 
detection of sky structure. For this reason we interpret the measured power spectra from these images as upper limits to the polarized 
sky structure. Artifacts from one known bright source 3C200 are faintly visible as concentric rings centered on the source location. 3C200 
is in the upper center of the image just outside the main beam. 



Q or U to I leakage in EoR observations be controlled to 
< 10"^ 

Our observed upper limit to polarized sky structure in 
frequency channels of width 1 MHz is in the range of 1-3 
K at 300 < I < 1000, with weaker upper limits at smaller 
scales. This is shown in Figure [J] For the field we studied 
the polarized sky structure is an order of magnitude smaller 
than one would have predicted by scaling from the published 
350 MHz observations. 

The results reported in Ide Bruvn et al.l (|2006l l are at 
different locations on the sky, so a direct c omparison is not 
straightforward. The Jde Bruvn et al.ll2006l ) field is centered 
at 13hl4m +45 (J2000), while our field is at 08h26m +26. 
However, one way to attempt a comparison is to look at the 
WMA P measurement of p olarization emission on the two 
fields (|Hinshaw et alll2008l ). 

Smoothing the WMAP data to 3.7 degree resolution, 
the 13h field has a polarized temperature of 13/iK, while 
our 08h field is 17/iK. The noise is comparable to the differ- 
ence, so we estimate our field has a similar net polarization 
at high frequency. I n the Haslam 408 MHz intensity maps 
ijHaslam et al.lll982r ). the 13h field has a temperature of 20K 
which is similar to the 08h field at 18K. Using the available 
data the fields seem quite similar. 

Ide Bruvn et al.l l|2006l ) explain the observed polariza- 
tion structures as due to a magnetic field that is smooth 
across their field, combined with an electron density that 
varies. This causes structure in the Faraday depth allowing 
the polarized sky structure to exceed the total intensity vari- 



ations. We see no reason to think our field would be different 

in t his regard. 

IVinyaikin et al.l l|l996l ) provide a model of diffuse Galac- 
tic polarization toward the North Celestial Pole (NCP). 
They find polarization \/Q2 + m of 3.5 ± 1.0 K at 88 
MHz (with a 24° FWHM beam) and 2.15 ± 0.25 K at 200 
MHz (8° FWHM)Q Their result is difficult to compare to 
ours, for several reasons. First, there is no overlap in {u,v) 
coverage between their single-dish (or phased antenna ar- 
ray) measurements and our interferometer results (which 
have no zero- length baselines). Second, one must interpo- 
late between fre q uenci es to compare to our 150 MHz results; 
IVinvaikin erall l|l996h can fit physical models with polar- 
ization nulls between 88 and 200 MHz, but the number and 
location of such nulls is not uniquely determined. Finally the 
observations are in different parts of the sky, and one would 
not expect geometry-dependent properties such as the fre- 
quencies of polarization nulls to be exactly the same. 



10.1 Depolarization Models 

Several depolarizing effects have been predicted at low 
frequencv(.Buru.l966|; ISokoloff et a"l]|l998l : IVinvaikinll2004l ). 
These include beam depolarization due to field gradients. 



^ Note that the numbers in Table 1 of lVinvaikin et al] lll996h are 
actually 2\/ + due to their normalization convention for 
polarized brightness. 
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Figure 6. The 3D polarization power spectra based on our ob- 
served visibilities. These are estimated from the power of maps at 
fixed rotation measures. The transverse scale has been converted 
from angular to comoving units at z = 8.5. The parallel scale is 
converted from ip according to Eq. IIA38I I. The contours are of Ap 
(see text), and start at 0.0069 K'^ in the left middle near fc ~ 0, 
increasing logarithmically by 10 per contour level. On large scales, 
fcx < O.S/iMpc"^ our observations set an upper limit equivalent 
to a few Kelvins per In k in 3D space. 

or sub-beam structure, and depth depolarization due to co- 
location of the rotating and emitting regions. We now con- 
sider each effect in turn as possible explanations for the low 
level of polarization at 150 MHz. 

First we conside r beam depolarization (|Burnl Il96d ). 
Ide Bruvn et all (|2006t ) reported RM gradients of VRM=2-4 
radm-^deg-^ At 150 MHz (A = 2 m) a gradient this size 
would result in a periodic modulation of the polarization 
vector field at 

I = VRM ~ 1000 - 2000. (6) 

TT 

We could detect such a periodic variation, but do not see 
such a feature. If the spatial gradients were substantially 
larger in our field, the periodic structures could occur on 
scales smaller than we constrain. However, this would re- 
quire a RM gradient larger than 20 radm"'^ deg~^. Beam 
depolarization due to presence of a smooth RM gradient 
seems an unlikely explanation for the low polarization. 

Alternatively, beam depolarization could arise if there 
was a small scale random magnetic field. An equipartition 
magnetic field in the galaxy could be several /xG, but this 
is a factor of several larger than the coherent field. The co- 
herent field is observed to contribute an RM=5.9 radm^^ 
to the pulsar, which results in a 23 radian rotation of the 
polarizati on vector. The dispers i on measure DM= 19.4751 
cm~ pc (IManchester et al. I l2005l : lManchesteilll974 ). results 
in a mean electron-weighted magnetic field along the line of 
sight field of (By) = 1.232RM/DM = 0.37^G. This coherent 



part of the field is lower than the equipartition value making 
plausible the idea that there is fine structure in the field. 

Substantial beam depolarization will occur if the small 
scale field structure produces variation of rotation of order 
unity across the synthesized beam. Given the current sensi- 
tivity at GMRT, structures at I ^ 10^ are unconstrained. 
A statistically isotropic field on scales smaller than that 
would have a correlation length 10* shorter than the line 
of sight. Random rotation measures add stochastically, so 
the net angle would be equivalent to a 100 times weaker, 
but coherent field. So a random magnetic field which has 
five times the field strength of the coherent field and a cor- 
relation length 10* times shorter than the rotation screen 
would lead to be am depolarization. Tak ing an electron scale 
height of 1 kpc l|Cordes fc Lazid |2002| ). which at b = 30° 
is 2kpc in projection, would correspond to a fine-structure 
scale of ~ 0.2pc. A magnetic field of 3 fiG on these scales 
would suffice to depolarize, which is in the plausible range of 
galactic magnetic fields. The spatial scales are a bit smaller 
than the outer scales o f 1 pc proposed in turbulence models 
tHaverkorn et al1l2008l ). but may still be plausible given the 
general nature of these arguments. At our galactic latitudes, 
the outer scale may also be substantially larger. 

We now take the analysis of fine scale beam depolariza- 
tion one step further by assuming a turbulence power spec- 
trum. We consider a Faraday thin emitting region with a ro- 
tating medium in the foreground. We model the RM with a 
Kolmogorov-like p ower spectrum on small angular scales, i.e. 
ARM^(6I) oc e^^^ (|Haverkorn et al.ll2008l ). In such a model 
the polarization direction is coherent over patches of size 6^ 
where ARM(6lc)A^ ~ 1. For the Kolmo gorov spectrum this 
implies 6c oc X~^^^^. The peak of the observed polarized 
power spectrum then s cales as ^peak ~ 6^^ oc A^^^®. Ob- 
servations at 350 MHz l|de Bruvn et al.ll2006l ) typically find 
polarized emission at scales of order 10 arcmin. There do 
not appear to be published power spectra of these polariza- 
tion maps but if the Q and U Stokes parameters change sign 
on this scale, as o ne would expect from t he abundant depo- 
larization canals (|Haverkorn et all 12004 ). then these struc- 
tures contribute mainly at I ~ 10''. Scaling to 150 MHz 
would move this power up to smaller scales by a factor 
of (350/150)^^/^ ~ 8, or Zp^ak ~ 10*. This reduces the 
tension with our observations, which have upper limits on 
of a few tens of K on these scales (see Fig. |4]). 
This is close to the expected signal from scaling from higher 
frequencies; the additional presence of depth depolarization, 
or different turbulence parameters in our field would elimi- 
nate any tension. 

Next we consi der depth dep olarization, which was al- 
ready discussed by iBurnI l|l966l ) . This occurs because the 
emission and rotating media are mixed, and some rotation 
occurs within the source. This effect has been reported at 
1.4 GHz for observations looking right into the Galactic 
Plane, where rotation measures can exceed 1000 radm~^ 
(|Uvaniker et al.ll2003l ). Our regime has a similar product of 
RMA^, so it is tempting to draw analogies. For this part 
of the discussion we assume a uniform magnetic field. The 
synchrotron emission arises from relativistic electrons, which 
make up a small fraction of the total electron count. How- 
ever, the Faraday rotation is caused by all free electrons. 

The strength of the depolarization depends on the pro- 
file of the emitting region. An example of a rapidly depo- 
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larizing model is a Gaussian emission region embedded in a 
medium with uniform rie and B. This might be a reason- 
able model of an emission region due to cosmic rays diffus- 
ing across magnetic field lines. Some prominent structures 
in the sky, such as the Cygnus loop, may be examples of 
such sources. This results in a Faraday dispersion function 
that is a Gaussian 
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The observed polarization vector is 
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This model has an exponential suppression with argument 
oc A'', so a factor of ~ 2 suppression at 350 MHz would 
be a factor of 10*^ suppression at 150 MHz. Using the dis- 
persion and rotation measures to the pulsar, an emission 
structure at a distance similar to that of the pulsar, with a 
Gaussian width 10% of the extent of the rotating medium 
would produce an order of magnitude attenuation of the net 
polarization between 350 and 150 MHz. 

However it is unlikely the electron density follows the 
precise Gaussian profile needed to achieve the deep cancel- 
lation of the polarization from different regions. We now 
consider an alternative scenario, which would lead to much 
less depth depolarization. Assume the emission spatial struc- 
ture traces the ionized gas structure perfectly. Two scenarios 
could lead to depolarization: (i) the ionized and relativistic 
populations are uniformly distributed, but the magnetic field 
is confined to flux tubes; or (ii) the ionized and relativistic 
components trace each other (for example, if they were car- 
ried by the same magnetic field). In both cases, the Faraday 
dispersion function is a top-hat, i.e. F{ij)) = Pq/I^oI between 
and (jio. This leads to a net polarization 
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This leads to an oscillating depolarization with a ~ A~ en- 
velope, which at 150 MHz gives an additional factor of only 
5 suppression compared to the 350 MHz data, less than the 
increase due to the spectral index. In this model one would 
expect 5-10 K polarized emission, which is not observed at 
I < 10^ 

The emission region could also be only partially over- 
lapping with the rotation region. Many complex permuta- 
tions are plausible. The simple models presented above are 
meant to illustrate that strong depth depolarization is pos- 
sible, given the known parameters. 

10.2 Prospects 

To further study the nature of the depolarization, multi- 
wavelength polarized observations of the same field are 
needed. Repeated observations during daytime would give 
ionospheric rotation measure variations that could separate 
beam squint from intrinsic polarization. Much more sensi- 
tive polarized maps then could be made. Such observation 
would likely concentrate on fields that show substantial po- 
larized fiux. 

Several steps could also be implemented in the future. 
Just as parallactic angle rotation allows us to decompose the 
polarization components of the pulsar, one could use time 



variable ionospheric rotation measure to unleak the sources. 
The ionospheric Faraday rotation changes through the day- 
night cycle, as we see in Figure [T] After a few months, the 
same point on the sky will be viewed through a substantially 
thicker ionosphere, and the Faraday depth of all sources will 
be modulated by an ionospheric contribution A(/). The ap- 
parent angle of polarization will change, while unpolarized 
light stays unchanged. This allows a separation of polar- 
ized and unpolarized components which does not require 
modelling of beam squint. Specifically, in circularly polar- 
ized baselines, a two epoch observations at 1,2 measures 
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We observe the two visibilities Vi,2 at the same azimuth and 
elevation, so knowing allows us to solve for the sky and 
leakage visibilities. 

To test the scenario of small angle beam depolariza- 
tion, one could compare the rotation measure of pulsars in 
lobular clusters. A sea rch of the ATNF pulsar database 
Manchester et al ] |2005') did not turn up any close pairs. 
M15 and 47 Tuc have multiple pulsars separated by arc min- 
utes which would be predicted to have fractional variations 
in RM of order 10% in a beam depolarization scenario. 

As far as the search for high redshift reionization is 
concerned, one wants fields with low Galactic polarized flux. 
The upper limit to polarization we report indicates that suit- 
able fields are available. 

Other steps to reduce the polarization leakage could 
include mosaicking a wide field, where the angular patterns 
would cancel, as was done in the work of Ide Bruvn et al.l 
(j2006l ). In addition, a raster of observation of the pulsar 
would allow measurement of the polarization calibration into 
the skirts of the main beam. 

Given the low level of polarized sky structure we report, 
along with the estimates of polarization leakage at GMRT 
presented in Figure |3l we can estimate that even a modest 
polarization calibration will allow measurement of EoR sky 
structure. The one sigma upper limit to sky structure mea- 
sured as a 3-D power spectrum [k^ Pp{k) at scales 
relevant to EoR (fcx ~ fc|| ~ 0.1/iMpc~^) is around 2 K, 
and the I to Q, U leakage at the edge of the main beam is 
0.03. Assuming that the Q, U to I leakage is also around 
0.03 the polarized sky structure on this field could create 
apparent sky structure that would interfere with detection 
of the EoR signal at the 60 mK in 3-D. This is the level 
of polarization artifact that may occur at the edge of the 
main beam. Leakage is smaller near the center. A factor of 
10 primary beam leakage correction would be sufficient to 
bring potential polarized leakage levels below the < 10 mK 
expected EoR signal. Polarization calibration made by ras- 
tering the pulsar across the main beam, along with mosaic 
observations may allow the beam squint terms to be mea- 
sured and corrected for to achieve this additional factor of 
ten reduction in the polarized emission artifacts, allowing a 
measurement of large scale structure during the Epoch of 
Reionization. The required calibration should be attainable 
using the pulsar technique. There are caveats however; for 
example, if there is much stronger polarization in neighbor- 
ing fields then polarization leakage through the sidelobes or 
diffraction spikes may be an issue. As we have seen in this 
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paper, bright unpolarizod sources outside the beam can ap- 
pear polarized, but at zero Faraday depth. 



11 CONCLUSIONS 

The GMRT EoR data has resulted in the strongest upper 
limits to polarized emission at 150 MHz to date. The ob- 
served polaxization is an order of magnitude smaller than 
had been expected from higher frequency observations. We 
note that comparison was not made on the same fields. We 
estimated the expected beam and depth polarization effects, 
but these are model dependent and poorly constrained by 
the data. Within the uncertainties, either mechanism could 
account for the lack of observed polarization. If beam de- 
polarization were the cause, it would imply a significant 
small scale magnetic field with structure at / > 10, 000. Oth- 
ers have reported depth depolarization at similar RMA'^, so 
depth depolarization may play a role here as well. When 
compared to the sensitivities required for EoR observations, 
the upper bounds we report are still two orders of magni- 
tude larger than the sensitivities required. The measured 
polarization leakage at GMRT is less than 10%, so a mod- 
est leakage correction, accurate to 10%, should make EoR 
observations feasible. The low polarized emission we report 
is promising for the search for Reionization. 
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APPENDIX A: POWER SPECTRUM 
CONVENTIONS 

We report our final results in this paper in terms of (lim- 
its on) polarized power spectra. This is appropriate because 
the theoretical predictions for 21 cm fluctuations from the 
reionization epoch are usually presented as power spectra, 
either 2-dimensional (C;) or 3-dimensional [J'(fe)]. An addi- 
tional advantage is that unlike many other statistical mear 
sures, power spectra can be measured unambiguously over 
some range of I (or k) from data that are missing zero-length 
baselines or contain significant instrument noise. The pur- 
pose of this appendix is to summarize the definition of the 
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power spectrum for both total intensity and polarization, 
and show how they can be derived from visibilities. 



Al 2D power spectra: total intensity 

We first consider the 2-dimensional case, e.g. where we con- 
sider diffuse emission at only one frequency. The intensity 
I{n) is then a function of angular position n. It can be 
Fourier-decomposed as 

I jH^im'--. (Al) 



Here I is the Fourier wavenumber traditionally used in cos- 
mology; it is related to u = {u, v) (which has units of waves 
per radian) hy I = 2nu. Finally I is the mean intensity of 
the field. 

The power spectrum Ci of a diffuse field is a statistical 
measure defined by 

{r{l)iil')} = {2nf Ci5^''Hl - I'), (A2) 

where 5'^' is the 2-dimensional Dirac delta function. It is 
related to the correlation function of the field by use of 
Eq. (IXTll . 



{[I{n)-I][I{n')-I]} = 



(2^) 



:Cie 



(A3) 



In particular, the RMS intensity fluctuations in the map can 
be written as 



■C, 



l^Ci dl 
~2^T' 



(A4) 



where in the last equality we have integrated over the direc- 
tion of I and noted that the area element d?l in the Z-plane 
is 2-Kldl. This integral is the reason why I^Ci/{2-k) is a com- 
monly used quantity: it describes the contribution to the 
intensity variance coming from a particular range of Fourier 
modes (specifically, per In/). 

The power spectrum can also describe the RMS fluc- 
tuations of a sky map smoothed by a beam. For example, 
convolving /(n) with a Gaussian beam of FWHM 9^ to gen- 
erate a smoothed map 7sm(Ti) is equivalent to multiplying its 
Fourier transform by the Fourier transform of a Gaussian: 



/sm(0 = e 



-S^I''/(161n2) 



hi)- 



(A5) 



The power spectrum of the smoothed map is then 
jg-*b' /(I6in2)j2 ^jjj^gg ^jjg pQ-^er spectrum of the original 
map, so the RMS fluctuation of the smoothed map is 



{[h^{n)-lY)= I ^e- 



e^;^/{8in2) di 
I ' 



(A6) 



We finally relate the power spectrum to the variance of 
visibilities observed on an interferometer. We first consider 
the visibility V{u) observed on an interferometer with base- 
line u (separation transverse to line of sight in units of the 
wavelength A), primary beam FWHM 6b. For our angular 
coordinate n, we take the origin to be at the center of the 
primary beam. A point source at position n with flux F (in 
Jy) will give a visibility 



V{u) 



-(41n2)|n|'^/9i: 



(A7) 



where the first factor gives the fringe phase and the second 



is the primary beam suppression. In a diffuse field, the ob- 
served visibility will be the integral of the intensity over the 
field of view. 



V{u) 



'J(n)d n. 



(A8) 



Using Eq. \A1\ . and noting that the rapid modulation of 
the fringe phase implies that the mean intensity / gives no 
contribution, we find 



V{u) = 



dH , 



m 



-2-Kiu-n -(41n2)|Ti|^/9J iln -.2 



d n 



(2,)» 

The variance of this visibility can be found from Eq. (IA2 



-27ru|^/(161n2) 



(A9) 



{\V{u)f) = (^)'/ d^zce-^'l'-^^"!''''*'"^'. (AIO) 
Thus the variance of the visibilities is directly related to 



the power spectrum at values of I within 



of 2tvu. 



In the practical case that the fringe spacing is much smaller 
than the primary beam (Sb|ii| 3> 1) and the power spectrum 
varies slowly with scale, we may approximate Ci by its value 
at 2nu and write 



8 In 2 



Ci—2ttii 



(All) 



This implies that the power spectrum can be determined in 
terms of the visibilities: 



F 

— Ci\i=2tvu ~ 
ZTT 



16 In 2 



\V{u)\'}. 



(A12) 



In our case, if the visibilities are calibrated in Jy, then 
the intensity /(/) has units of Jysr~^ and the power spec- 
trum is determined in Jy^ sr~^. However we want the power 
spectrum in units of K^, which means that the visibilities 
have to be converted to Ksr units. The conversion factor is 

lJy = 0.00145 (i^)\sr. (A13) 

Inserting this conversion factor into V{u) in Eq. (|A12|l yields 

— ( 

2tt 



I r ! _ooc,,, in-sl^r / 150 MHz 

— O; k2 \i=2ttti — 2.66 X iU —J- 

(A14) 

If we refer the beam to the intrinsic FWHM of the GMRT 
primary beam, 3.3° or 0.058 radians, this becomes 



27r 



Ci |i = 2Trii 



u\ 3.3° \ / 150 MHz 



11.9 



\Vjy{u)\'}. 

(A15) 

Unlike e.g. the microwave background or extragalac- 
tic large scale structure. Galactic foregrounds are not sta- 
tistically isotropic. Nevertheless, one may still consider 
the power spectrum in a particular region of sky using 
Eq. (ETill . 



A2 2D power spectra: polarization 

The same Fourier decomposition used for the total intensity 
can also be applied to linear polarization. One may decom- 
pose the linear polarization in analogy to Eq. (|A1|I as 



Q(n) = 



dH 

{2^Y 



(A16) 
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and similarly for U. 

In studies of the microwave background, it is common 
to decompose the polarization into so-called i?-mode and B- 
mode components, in which one rotates the reference axes 
of the Stokes parameters to be aligned with the Fourier 
wavevector I: 



Q{1) 
U{1) 



cos{2ai)E{l) 
sm{2ai)E{l) - 



-sm{2ai)B{l), 
cos{2ai)B{l), 



(A17) 



where ai = tan ^ ( ly/lx. 
IZaldarriaga fc SeliakI Il997l : 



1996'; 



Stebbins '1996'; 'Soliak' ^199^ 
Kam ionkowksi ct al. 1997)0 
One can then define power spectra Cf and in analogy 
to Eq. HA2p . This is useful because in the case of the mi- 
crowave background the source of polarization is Thomson 
scattering, and hence the polarization direction is correlated 
with spatial morphology, leading to different power spectra 
in the E and B modes. For low- frequency Galactic emis- 
sion, however, Faraday rotation makes it unlikely that the 
observed polarization angles would have anything to do with 
the spatial morphology. Therefore we instead work in terms 
of the sum of the power spectra, Cf = Cf + CP , which 
incorporates both Stokes parameters. Then we have 

{Q*{l)Q{l') + U'{l)U{l')) = (27r)'cr5'''(« - I'), (A18) 

and the variance of the polarization, analogous to Eq. (IA4|) . 
is 

{Q\n) + U\n))^ /^f- (^19) 

Thus PCf /{2Tr) can be viewed as the contribution to the 
polarization variance per In/. 

The same methodology for estimating the intensity 
power spectrum from the variance of the visibilities can be 
applied to the polarization. We find: 



lul 3.3° \V 150 MHz ^ " 



11.9 



{\VLnMur). 

(A20) 



A3 3D power spectra: total intensity 

The predicted 21 cm signal is 3-dimensional in the sense 
of having structure in both the angular and radial (fre- 
quency) directions. The statistic most commonly computed 
by theorists is power spectrum, P{k,fj.), and the power 
per logarithmic range in wavenumber, fc^P (fc,^i)/(27r^) 
l|Bharadwai fc Alill2004l : iBarkana fc Loebll2005l ). To define 
this statistic, we begin by decomposing the intensity into its 
Fourier modes in analogy to Eq. (|A2Hl : 



J(r) 



{2nr 



7(fe)e'' 



(A21) 



here r is a 3-dimensional comoving position and k is the 
wave vector. Then a power spectrum is defined via: 

{i*{k)I{k')) = {2TTfP{k)S^^'>{k - fc'). (A22) 



Even in the case of a cosmological signal, the power 
spectrum may depend on the direction of fc relative to the 
line-of-sight vector which we take to be 63. This is because 
the observed frequency of emission depends not just on the 
location of the emitting hydrogen gas, but also on its radial 
velocitylfl We therefore introduce 



cos Oh ~ -r- 



k ■ 63 



(A23) 



which is between 
fluctuations is 



-1 and -1-1. Then the total variance of 



Ak Ank^P{k,ii) 



k 



2 27r2 ' 
(A24) 

where we have used that the volume element d^fe is 
27rfc^d/idfc after integrating out the azimuthal angle. For 
an isotropic power spectrum, P{k,^) — P{k), we can thus 
interpret k^ P{k)/{2-K^) as the variance in intensity per loga- 
rithmic range in wavenumber. Often one denotes A'^(fc, /i) = 
k^P{k, fi)/ {2-K^). The power spectrum is also sometimes 
written in components transverse to and radial along the 
line of sight, P(A;j_,fc||). 

In order to connect P{k, jj.) to observations, we need 
to find the conversion from the comoving coordinates r to 
the observables, angular position n and frequency v. In the 
direction transverse to the line of sight, the relation is rx = 
Dn, where D is the comoving angular diameter distance to 
redshift z — A/ (2 1cm) — 1. In the direction along the line of 
sight, the relation is more complicated: 



dry dr||/dz 

'di7 ^ AvjAz ^ -vl(\ + z) 



(cH)- 



l + z 
cHv ' 



(A25) 



where c is the speed of light, H is the Hubble constant at red- 
shift z, and V is the observed frequency. We can then define 
the radial conversion factor C = (1 -I- z)/{cHv). It follows 
that the comoving coordinates are related to observed co- 
ordinates via (rx,?"!!) <-> [Dn,—C{y — vo)]i where vq is a 
reference frequency (the zero point of is arbitrary). 

We can now construct the 2-dimensional power spec- 
trum, which can be measured via Eq. (|A12|I . in terms 
of the 3-dimensional power spectrum. We construct a 2- 
dimensional image l2o(n) with bandwidth Av, 



I{n, Au 



rCAv/2 



— - / /(Z?n,r||)dr|| 

<-'^J^ J -CAv/2 



(A26) 



Using Eq. (|Aip . we can find the Fourier transform of the 



^ Some early papers denoted the rotated polarizations by G and 
C, or ® and instead of E and B, or used different normal- 
izations. Also note that this decomposition of polarization also ^ An analog ous situation occurs in galaxy redshift surveys 
applies to other spin 2 fields such as gravitational lensing shear. llKaiseilll987t) . 
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2- dimensional map. For i 7^ 0, 
/2d(Z) = J d^n I{n)e'" " 



CAv 

1 



CAix/2 



-CAv/2 
CAv/2 

-CAv/2 



dr|| I[Dn, r||)e '' "^ 



drii 



d^fe I dfcii ^, 



(27r) 



L/(fc) 



2^(2)^ 



x(27r)"<5^^^(Dfc 



dm 



d^fcidfci, ~ 



(27r)3 
i)e"=iri! 



I{k) 



CAi/ku 



(A27) 



where sincx = s'mx/x. In the last equality, the d^k± in- 
tegral is trivial due to the (5-function, and the ry integral 
gives the sine function. By taking the variance we extract 
the 2-dimensional power spectrum, 

,n 1 f dkn fl \ yCAiykii 
Q^° = ^/^P(^,fc„)sinc^^. (A28) 

Note that the amount of power in the 2-dimensional map 
depends on the bandwidth Au. This is a direct consequence 
of the 3-dimensional nature of the EoR signal: fluctuations 
with |fc||| > {CAu)~^ are cancelled by the combination of 
frequencies. 

It is possible to measure the full 3-dimensional power 
spectrum by stacking maps at several neighboring frequen- 
cies with a phase shift to avoid cancellation of frequencies, 
i.e. 



hnin, t) = 



Au 



e '^"^"/(n, u) du, 



(A29) 



where r is the phase shift in units of cycles per MHz. This 
is now a complex map but this causes no difficulty. The 
power spectrum of hoirijT) can be obtained by the same 
methodology as Eq. (|A28p . yielding 

(")=o^y ^^U-'^v^'^ — 2 — • 

(A30) 

Thus this stacked map is sensitive to fluctuations with ra- 
dial structure and wavenumbers in the range fc|| ~ 2nT/C, 
witihin a width (5fc|| ~ {CAu)~^. If the power spectrum is 
slowly varying as a function of fcy over this range, one may 
use the fact that J sinc^a;da; = tt and approximate 



crir) 



1 pfl ^\ 

D^CAv \D' C 



(A31) 



By combining with Eq. (IA14|I . this allows one to estimate 
the 3-dimensional power spectrum from frequency-stacked 
2-dimensional visibilities. 



A4 3D power spectra: polarization 

The machinery used to describe intensity fluctuations can 
also be used to describe polarization fluctuations. Polarized 
emission from EoR is expected to be intrinsically very small 
and if present would be strongl y affected by Ga lactic Fara- 
day rotation (although see ■Babich fc LoebH2005l ). Polarized 



emission from the Galaxy is a potentially severe foreground, 
and we characterize it by describing its power spectrum us- 
ing the same conventions used for the EoR signal. To do 
this, we decompose the polarized flux density into Fourier 
modes just as was done for the intensity: 



Q{r) = 



(2^ 



(A32) 



and similarly for U . One may define a polarized power spec- 
trum via 



{Q*{k)U{k')) = (27r)^PQc(fc)5(^'(fe - fe') 



(A33) 



and similarly for Pqq, Puq, and Puu- Since Q and U are 
real it follows that Pqq and Puu are real and non-negative 
and satisfy Pggik) = Pqq(— fe) and Puu{k) = Puui—k); 
also Pguik) ~ PuQ^—k). A decomposition of the polar- 
ization tensor into E and B modes is possible, just as in 
the two-dimensional case, but is probably not interesting 
if the radiation has propagated through a Faraday thick 
medium. We therefore focus on measuring the power spec- 
trum Pp{k) = Pqq + Puu- 

{Q'{k)Q{k') + U''{k)U{k')) = {2nfPp{k)S^^Hk ~ fc')- 

(A34) 

The 3-dimensional polarized power spectrum can be 
obtained by stacking in analogy to Eq. (|A29|I . We define 
stacked 2-dimensional polarization maps. 



P2D{n,T) 



1 

a17 



I/0+AI//2 



1^0— A1//2 



'[Q(n, 1/) + iU{n, z^)] du. 

(A35) 



Repeating the steps that lead to Eq. (|A3ip gives, in the 
polarization case, that the 2-dimensional power spectra of 
the complex maps P2D(n,r) are 

1 



D'^CAv 



I 2nr\ 



2^Pqu 



I 2Tvr\ 



(A36) 



where 9 denotes the imaginary part. Note that for a signal 
that respects statistical parity invariance, the second term 
here must vanish, but this may not be true for the Galactic 
emission in a particular patch of sky. This term may be 
removed by averaging signals from positive and negative r. 

In the case of the cosmological signal in the unpolarized 
intensity, we identified fcy with the line-of-sight wavenumber 
of a mode in large scale structure. In the case of polarized 
Galactic synchrotron radiation, the rapid variation of the 
signal with frequency is instead due to Faraday rotation. 
In fact, Eq. HA35|I is (at least over a limited range of fre- 
quency) an implementation of RM synthesis, with the ex- 



ponent e 



-2i<t>\^ 



replaced by e 



That is, it is the RM 



synthesis map at Faraday depth 



dv 
-TT— 7— — r 



d(A2 



2A2 



(A37) 



Using the conversions from r to fcy (fcy = 2-kt /C), this be- 
comes 



Cv 



l + z 



fell 



(A38) 



^ 4A2 " -icHX^ 
We thus conclude that the mathematical procedure to 
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extract the EoR power spectrum from a foreground-cleaned 
data cube is (in principle!) similar to that of computing a 2- 
dimensional power spectrum on RM synthesis maps, except 
that it is applied to the real total intensity / [via Eq. HA29|) ] 
rather than complex polarization Q + ill. In particular, a 
simple type of polarization-to-intensity leakage such as 

7obs(M,«,!^) = I{u,v,u)+'yQQ{u,v,u)+'yuU{u,v,u) (A39) 

will leak polarized structure at Faraday depth into the in- 
tensity map at wavenumber fc|| given by Eq. (|A38|I if (7q , "fu) 
is frequency-independent (or varies slowly with frequency), 
since it results in an apparent Stokes / signal with a si- 
nusoidal variation in v. In this case, the observed power 
spectrum would be 

PobB(fc) = P{k)+'y^QPQQ{k)+'ylPuu{k) + 2'yQ'yu^PQu{k). 

(A40) 

More complex types of frequency-dependent polarization 
leakage can occur in real EoR observations, especially if ob- 
servations at the same {u, v) but different frequencies are 
taken non-simultaneously [as is necessary at GMRT, since 
a given point in the (m, v) plane corresponds to a different 
physical baseline length at a different frequency]. One ex- 
ample would be polarized sidelobes that rotate relative to 
the sky due to the alt-azimuth mount. Another would be 
changes in ionospheric Faraday rotation between the epochs 
of observation. Finally the beam squint will vary across the 
bandpass. These frequency-dependent leakages are a concern 
because polarization at any Faraday depth can be injected 
into the Stokes / map at wavenumbers that do not obey 
Eq. dHHl. 
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